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BEHAVIOUR OF TALL BUILDINGS DURING 
THE CARACAS EARTHQUAKE OF 1967 (1) 


RESUMO 


Eistuda-se o comportamento de s5G edifícios e outros 
tipos de estrutura durante o sismo de Caracas de 1967. 
Calculam-se as frequências próprias e os coeficientes 
resistentes para forças horizontais e comparam-se com 
os danos observados. Em grande número de casos me- 
diu-se no local a frequência própria. Estima-se a inten- 
sidade do sismo em diferentes sonas a partir dos deslo- 
camentos máximos sofridos pelas estruturas e das forças 
necessárias para produzir os danos observados. Refere-se 
brevemente a influência do solo, Os resultados obtidos 
neste trabalho prestam valiosa contribuição para o julga- 
mento das regras de dimensionamento sísmico de edifi- 
cios altos de betão armado e em particular dos coefi- 
cientes sismicos actualmente adoptados. 


1. INTRODUCTION 


The study of the structural behaviour of buil- 
dings under strong earthquakes yields results of 
paramount importance for judging design rules. 
This is particularly true in the case of the Caracas 
earthquake of 29th July 1967, which struck a 
big town with more than 1000 buildings over 
10 stories. Due to the earthquake several buil- 
dings (5 of which more than 10 stories high) 


by  J. FERRY BORGES *, 
J. GRASES ** 
and A. RAVARA *** 


SYNOPSIS 


The behaviour o/ buildings and other structures 
during the 1967 Caracas earthquake is analysed. Static 
and dynamic characteristics, such as coefficients of 
resistance and natural frequencies, are computed and 
compared with the observed damages. In several cases 
natural frequencies were also measured on the site. 

The intensity of the ecarthquake in different areas is 
estimated on the basis of the maximum displacements 
suffered by the structures, the forces necessary to produce 
rupture and the amount of damage observed, The in- 
Juence of soil conditions is briefly discussed, 

The results obtained may contribute in assessing the 
adequacy of the design rules concerning reinforced con- 
crete tall buildings and, particularly, of the seismic 
factors usually adopted in the codes. 


collapsed and many others were considerably 
damaged. 

In the survey prepared by the Ministerio de 
Obras Públicas [1], 2300 cases of damage were 
recorded in the Caracas area, from which 27º/, 
were considered important. One half of these 
cases concern one or two story houses and the 
other half higher buildings. Sixty-eight struc- 
tures, 17 of which were buildings, collapsed or 
were damaged beyond repair. 


(1) Paper presented to the IV World Conference on Earthquake Engineering, Santiago de Chile, Jan. 1969. 
* Associate Director, Laboratório Nacional de Engenharia Civil, Lisbon. 
“* Professor, Instituto de Materiales y Modelos Estructuráles, Facultad de Ingenieria, Universidad Central de 


Venezuela, Caracas. 


*** Head, Applied Dynamics Division, Laboratório Nacional de Engenharia Civil, Lisbon. 
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In the present report 51 cases are referred. 
A quantitative study of the structural behaviour 
is presented for 22 structures, 14 of which are 
buildings (9 having 10 stories or more). Four of 
the collapsed buildings are studied. 

Each structure is analysed on the basis of 
information obtained in the direct inspection 
of the site, on results of tests of materials per- 
formed at the IMME [2] (Instituto de Materiales 
y Modelos Estructurales) of the Universidad 
Central de Venezuela and on the available plans. 
It was possible in each case to compute the 
values of the horizontal forces which if applied 
at a given level and in a given direction would 
produce collapse. By dividing these forces by the 
weight of the structure above the considered 
level, coefficients of resistance were obtained. 
These resistance coefficients can be directly com- 
pared with the seismic coefficients prescribed in 
the codes if they are computed according to the 
same design rules. 

Also in each case the natural frequencies of 
each building are computed and compared with 
experimental values when available. 

The intensity of the earthquake at the different 
locations (expressed in power spectral density of 
acceleration) is estimated on the basis of the 
maximum displacements undergone by the struc- 
tures, on the basis of the ultimate forces they 
can resist and by judging the observed damages. 
The increase in seismic intensity due to the in- 
fluence of soil conditions is also briefly studied. 

By comparing the observed damages with the 
results of the computations, conclusions are 
drawn concerning the intensity of seismic forces 
and the adequacy of design rules. 

No general description of the earthquake and 
of its consequences is attempted in the present 
paper. Particular attention is called to recent 
publications on the matter [3 to 10). 


2. BASIC ASSUMPTIONS 


The interpretation of the results concerning 
structural behaviour is based principally on the 
comparison between the amount of observed 
damage and the coefficients of resistance to hori- 
zontal forces. The natural frequencies correspon- 
ding to the type of structural behaviour assumed 
for computing the coefficients of resistance are 
also taken into consideration. 


Damage is classified into 4 groups: collapse, 
important damage, small damage, no damage. 
This classification tries not only to express the 
amount of damage observed but also to express 
a subjective probability of collapse in correspon- 
dence with the observed damage. The structures 
with no damage and small damage are conside- 
red to have a small probability of collapse. Those 
structures assumed to have a probability of 
collapse near to 1 are classified in the important 
damage group. Thus the structures with impor- 
tant damages present a clear mechanism of col- 
lapse and local damage or residual displacements 
showing that their ductility is practically 
exhausted. 

For determining the static horizontal forces 
necessary to produce the collapse of the reinforced 
concrete structures, and thus the resistance 
coefficients, the following hypotheses were adop- 
ted concerning the type of structural behaviour 
and the ultimate strength of reinforced concrete 
elements. 


2.1 — Structural behaviour 


The observation of the damaged buildings 
suggested the 6 main types of structural behaviour 
indicated in fig. 1. Thus the ultimate forces 
necessary to produce collapse were computed 
according to these idealizations. 


C) 


PLASTIC 
HINGES 


Í 


Fig. 1 — Tipes of structural behaviour 


Structural behaviour type 1 corresponds to 
the formation of plastic hinges at the top and 
at the bottom of the columns at one of the stories, 


According to ultimate design principles the 
horizontal force which produces collapse is 
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obtained by adding the ultimate forces of the 
different columns. This type of structural behaviour 
occurs when the beams of the floors are conside- 
rably more rigid and resistant than the columns 
or when the walls at the stories near the considered 
one can prevent the rotation of the joint blocks. 
This type of behaviour often occurred in buildings 
with a small amount of partitions at one level 
(in general the ground floor). 

Structural behaviour type 2 corresponds to 
a collapse mechanism affecting two stories, with 
formation of plastic hinges in the columns and 
in the floor structure (beams or slabs). This type 
of behaviour occurs when the structure of one 
floor cannot resist the moments transmited by 
the columns. This is a very frequent case in 
Caracas, where structures with transverse beams 
only are very often used. Generalizing this type 
of behaviour it would be possible to consider 
the formation of plastic hinges in more than two 
stories but such schemes are no more realistic. 
In fact plastic hinges can only form at several 
stories under too large horizontal displacements. 

Structural behaviour type 3 corresponds to 
the elastic behaviour of framed structures. When 
the buildings are formed by parallel frames of 
unequal stiffness, accurate results are only obtai- 
ned if compatibility is imposed to the horizontal 
displacements of the different frames at the dif- 
ferent levels. This type of frame behaviour only 
occurs when the structural contribution of parti- 
tions is small. 

Structural behavior type 4 corresponds to the 
association of shear walls with a framed structure. 
This type of behaviour is studied using the 
usual frame analysis programs but introducing 
segments of infinite rigidity to reproduce the 
dimensions of the joint blocks. 

Structural behaviour oftypes 5 and 6 corres- 
pond to cantilevers with distributed mass and 
concentrated mass at the top, respectively. 

It must be noted that the type of structural 
behaviour may change as a function of the 
intensity of the seismic loading. 


2.2 — Ultimate strength of reinforced concrete 
elements 


The ultimate strength of reinforced concrete 


elements under pure bending or combined axial 
compression and bending is computed according 
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to the Recommendations of European Committee 
for Concrete [11]. These Recommendations were 
adopted as a basis for the Portuguese Code for 
Reinforced Concrete Structures [12] and thus 
the tables prepared for this code were used. 

The basic assumptions consist in considering 
for concrete a stress-strain diagram defined by 
a parabola with the vertex at point (7, , 2º/s0) 

7 
with 0) = E 7,» being the characteristic (5 º/, 
fractile) of the 28 days cylinder strength of 
concrete. This cylinder strength is taken as 
85º» Of the cube strength. 

For mild steel the elasto-plastic stress-strain 
diagram is defined by a yelding stress 0“ = 
= e , €;, being the characteristic yielding 
stress of the steel. For work-hardened steel the 
diagram is curved and defined by the character- 
istic limit of proportionality at 0.2º/,, divided 
by the minoration factor of 1.15. The ultimate 
strain for steel is 10 º/,,. 

For simplifying the analysis of the data the 
concrete qualities are classified in 4 categories: 
B180, B225, B250 B300, the reference numbers 
indicating the characteristic cube strength. The 
steel qualities are designated by A24 and A40 
corresponding to the characteristic yielding stress 
of 2400 kg cm? and the characteristic limit of 
proportionality at 0.2º/, of 4000 kg 'cm?. 

The interaction diagrams of fig. 2 indicate the 
ultimate values of axial force, N, and bending 
moment M, for symetrically reinforced rectangu- 
lar columns with total percentages of steel, 
w,==0,5,1.0 and 2.0º/, for the following combi- 
nations of the mechanical properties of the mate- 
rials: i) A24, B180 and ii) A40, B300 (curves 1). 

The interaction diagrams represented by curves 

2 correspond to the ACI 318-63 Building Code 
(Chaper 19, Combined Axial Compression and 
Bending — Ultimate Strength Design). Conside- 
ring elastic behaviour and working stress design, 
curves 3 are obtained. Curves 3 refers to the 
allowable stresses currently used in seismic 
design, which are 33º/o higher that that those 
for permanent and live loads. 

A comparison of curves 1 and 2 shows that 
the CEB or ACI criteria are not very different. 
ACI criteaia are about 20 to 30º/, safer than 
the Portuguese Code of 1967. The figure empha- 
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sizes the inadequacy of working stress design 
for compined axial force and bending moment. 
In fact the ratio of the bending moments indi- 
cated by curves 1 and 3 vary between about 
1 and 2. 
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The safety concepts in which these combina- 
tions are based are not clear. From a probabilistic 
point of view it does not seem justifiable to 
consider earthquake loads multiplied by 1.4 in 
expression i) and by 1.25 in expression ii). 


M 1- CEB (PC.-1967) 
N 
| A 24 Ga=2400. 2100 kg/cm2 
a A 40 Ga +2 900 =3500kg/cm?2 


B 180 Gp SBSAEO. 102 kg/cm? 


B300 Ob= ASaaDO 170 kg/cm? 


* n4 e 
A 24 ty =2400kg/cm? 
DA 
A4O ty =4000kg/cm? 
A 40 R 
B 300 B1BO f. =085x180=153 kg/m 
B300 f. =085x300=255kg/m? 
(1 
3,- ELASTIC DESIGN, 
WORKING STRESSES 
| Vo A24 Ga=-2x1400=1870kgiem? 
niid WD ,=05% 3 
-— = 10% A4O Ga: x2400:3200kgicm?2 
— ) ,=2 0% 3 
B180 Gp=Lx 70-93 kgicm? 
50 60 3 


B300 Op=-S-xn15 =154kg/cm? 


Fig. 2 — Interaction diagrams for reinforced concrete columns 


2.3 — Load factors 


For comparing the coefficients of resistance 
with the seismic coeficients indicated in the codes 
attention has to be paid to the values of the 
specified load factors. In fact several codes take, 
for earthquake loading, load factors equal to 1 
but other codes prescribe load combinations 
where earthquake loads are multiplied by factors 
different from 1. This is, for instance the case 
of the 1967 revision of the American code 
«Recommended Lateral Forces Requirements» 
[13]. This code indicates that load resisting 
frames must be designed for the following loads : 


) U=1.40(D+LA+AE) 
ii) U=09D+1.25E 


where D — permanent load, L — live load, E 
earthquake load. 
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Anyhow, when comparing coefficients of resistance 
with seismic coefficients specified in this American 
code, the indicated values of the load factor have 
to be considered. Also if working stresses are 
used (even if increased by 33'/,) the ultimate 
bending moments obtained are in general much 
smaller than those determined in ultimate strength 
design. Thus the values of coefficients of resis- 
tance determined in the present paper can only 
be directly compared with the seismic coefficients 
of the codes that use a load factor equal to 1 
and ultimate load design. 


2.4 — Natural frequencies 


The interpretation of the results is based on 
the natural frequencies of the structures, which 
were computed and compared with experimentally 
determined values, 
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For computing the natural frequencies struc- 
tural schemes analogous to those used for deter- 
mining the coefficients of resistance were adopted. 

In the case of structures of types 1 and 2 
deformability of one or two stories only was 
considered, the remaining being assumed to be 
rigid. The behaviour of the deformable stories is 
assumed to be perfectly elastic and thus linear 
one-degree-of-freedom oscillators are obtained. 
For structures of types 3 and 4 the natural 
frequencies are determined by the Rayleigh 
method and in some cases also by modal analysis. 
In these cases the effect of non structural 
walls is disregarded and masses are concentrated 
at floor levels. In cases 5 and 6 agiven amount 
of foundation deformability is assumed in order 
to obtain values of the natural frequencies in 
agreement with the measured ones. 


3. STUDIED CASES 


The studied cases are summarized in Tables | 
to III, which refer to structures at Palos Grandes 
and Altamira in the town of Caracas (Table 1), 
at Caraballeda (Table II) and at Naiguatá and 
Playa Grande (Table III), fig. 3. At Palos Grandes, 
Altamira and Caraballeda there are deep layers 
of alluvial soil. At Naiguatá and Playa Grande 
the thickness of soft soils is supposed to be 
much smaller, 


nn 


| | | 10º 40' 
ros) enciren | | 


1050 


10º 40 


10º 30' 


10º" 25 


67"00' 66"50 


67º20' 67"10' 


Fig. 3 — Epicenter location. 


The qualities of the materials were experi- 
mentally determined in some cases [2]. When 
test results were not available the qualities 
indicated in construction plans were adopted. 
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67º00' 


When no indication was given in the plans the 
qualities A24 and B180 were assumed. 

The tables indicate the approximate orientation 
of the structures and the ratios H/L of height 
to longitudinal and transverse dimensions at the 
base. 

The amount of damage and the type of 
structural behaviour are classified as indicated 
in 2). The estimated importance of the structural 
contribution of the partitions is also included. 

The cofficients of resistance, and the computed 
and measured natural frequencies in each direction 
are finally indicated. Natural frequencies were 
measured with a Dr. Baule's B4h displacement 
meter produced by Hóttinger [14]. 

Fig. 4 to 13 show some of the collapsed or 
damaged structures. 

Fig. 14 compares the importance of damage, 
the coefficients of resistance and the computed 
natural frequencies at the Palos Grandes — Alta- 
mira area and at Caraballeda, Naiguatá and 
Playa Grande. 

Table IV summarizes the available data on 29 
cases that were only summarily observed and 
for which the coefficients of resistance could 
not be computed. The natural frequencies indicated 
for all these structures were measured in some 
cases before and in other cases after the structures 
had been repaired. Thetype of structural behaviour 
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indicated corresponds to the observed damage. 
Thus, when no damage was observed no indication 
is given. 

For the Residence Union and Mene Grande 
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(1) DETERMINED WITH REINFORCING SCAFFOLDING INSTALLED. 
(2) CORRESPONDS TO THE HORIZONTAL FORCE NECESSARY TO PRODUCE 
FAILURE IN THE DAMAGED COLUMNS. THERE IS AN IMPORTANT RESERVE OF 
RESISTANCE DUE TO THE SHEAR-WALLS AND COLUMNS NEARBY. 


TABLE |l - STUDIED STRUCTURES AT CARARALIEDA 
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(1) CORRESPONDS TO THE 7th STORY BUILDING REMAINING AFTER THE COLLAPSE 


OF THE UPPER STORIES. 
(2) REFERS TO THE qnd STORY. 
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TABLE IV - DATA ON OBSERVED STRUCTURES 


MEASURED TYPE OF 
. H/L NAT, FRE- STRUCTURAL OBSERVED DAMAGE 
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NO STRUCTURAL DAMAGE. SMALL, DAMAGE IN 


Ds esa = Se. 


ig. 5 — Mansion Charaima 


buildings (cases 41 and 47) which are H — shaped 
in plan, the torsion vibration mode was detected. 
This mode was particularly obvious in Residence 
Union. 

Fig. 15 relates the natural periods measured 
in longitudinal and transverse directions to the 
number of stories. 


4. INTERPRETATION OF THE RESULTS 


4,1 — Intensity of the earthquake 


For the interpretation of the results it would 
be of paramount importance to estimate the 


intensity of the earthquake in different areas 
and to study the influence of soil conditions. 
As no strong motion instrumental records are 
available this is a difficult task. 

A way to estimate the power spectral intensity 
of acceleration is based on the maximum displa- 
cements observed in some structures. According 
to previous results [15] the mean maximum 
displacement of a one-degree-of-freedom oscil- 
lator with 5º critical damping, expressed in 
cm, is approximately given by 

o My E. Vs 
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Fig 6 — Palace Corvin 
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Fig. 8 — San Bosco 


where f (Hz)* is the natural frequency of the 
structure, 5 (cm?s”!Hz) is the power spectral 
density of acceleration of the considered earth- 


(*) Hz — Herz — cyles per second. 
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Fig. g — Castillete 


quake and 700 cm*s *Hz is the intensity of 
the staddard earthquake (30 s of N-S El-Centro 
1940 record). Thus S can be estimated by 
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Fig 10 — Mene Grande 
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Fig. 12 — Macuto Sheraton Hotel 


The computations concerning the Macuto- 
Sheraton Hotel main building (Table II, ref. 
no. 14) indicate that for a horizontal transverse 
force equal to 0.22 of the weigth of the structure 
(force necessary to produce vyielding of the 
columns) the elastic displacement above the 4th 
story is of about 4.5 cm. Tests on reinforced 
concrete models to scale !5.; of the collapsed 
columns have showa that the amount of damage 
observed in the columns corresponds to residual 
displacement of about 3.5 cm (fig. 16 and 17). 
Thus the maximum displacement due to the 


3520 


“o Ma 


| ANA Raw 


Fig. 11 — Plaza Uno 


Fig. 13 — Pergola, Macuto Sheraton 


earthquake can be estimated as 8 cm. By using 
expression 2) and by considering the natural 
frequency of 1.4 Hz, a spectral density of 1400 
cm 25-*/Hz is obtained. This value is twice the 
indicated standard one. 

For the Macuto-Sheratton Pergolas (Table II, 
ref. no. 22 and 23) the indicated expression is no 
longer valid as the effect of gravity can no more 
be disregarded, as indicated by Husid [16]. By 
performing a dynamic non-linear analysis taking 
in consideration the effect of gravity, a maximum 
displacement of 24 cm is reached (for t=-13 5) 
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Kig. 14 — Relation between amount of damage, coefficient 
of resistance and computed natural frequency 


Fig. 16 — Model test of column 
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when Pergola 1 is subjected to the accelerogram 
of N-S El-Centro 1940 earthquake having its 
ordinates multiplied by V1.5. 

Another way to estimate the power spectral 
density of acceleration consists in computing the 
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Fig. 15 — Measured periods 


response of systems with a number of degrees 
of freedom sufficient to secure an accurate 
representation of the structural behaviour and 
to verify at what values of the spectral density 
of acceleration the observed amount of damage 
or collapse is attained. Computations of this 
type considering non-linear behaviour were per- 
formed for Mansion Charaima (Table II, ref. no. 
17) and show that for a power spectral density 


of 300 cmºs—*/Hz rupture of the columns is 
attained at the 7th story due to a whiplash 


effect, fig. 18. 

The dynamic analysis of Petunia II building 
(Table I, ref. no. 8) performed considering the 
linear behaviour of a 11 degree of freedom 


Fig. 17 — Macuto Sheraton column 
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system on an elastic foundation indicates that 
collapse by compression failure of the external 
columns occurs for a power spectral density of 
acceleration of about 500 cm? s”*/Hz. 

The seismoscope situated at Observatorio 
Cajigal (T == 0.7 s, 1 == 0,029) founded on rock, 
suffered a maximum acceleration of 200 cm*/s 


alluvial layers at Palos Grandes — Altamira with 
a thickness over 80m and of softer alluvial 
layers reaching to at least 70m in Caraballeda 
area. These thicknesses obtained through recent 
drilling [19] exceed the ones mentioned in pre- 
vious geological information [7]. 

Spectral analysis of microtremors registered 
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Fig. 18 — Dynamic analysis of Mansion Charaima 


[18]. This acceleration corresponds to a constant 
spectral density of about 150 cm?s-*/Hz over 
the frequency range 0-5 Hz. 

Finally comparisons of the coefficients of 
resistance with the observed damages, and the 
observed damages themselves give qualitative 
indications on the intensity of the earthquake. 

It can be concluded that the intensity of the 
earthquake was much amplified in the Cara- 
balleda area where an intensity considerably 
higher than the standard one was attained. This 
amplification is also apparent at Palos Grandes — 
Altamira but the intensity in this area is assu- 
med to be lower than 700 cm* s—*/Hz. In 
Caracas outside this area the intensity of the 
earthquake was even smaller. 

It must be noticed that, according to the 
magnitude of the earthquake (about 6.3 Richter) 
and the distance to the epicenter (70 km) [18], 
an intensity smaller then the one of El-Centro 
would be expected in Caracas. 


4.2 — Soil conditions 


Available data on soil conditions in the studied 
areas indicate the occurrence of consolidated 


by Grasses [20] at three different places of Palos 
Grandes — Altamira show a predominant fre- 
quency of about 0.7 Hz. This value, that is near 
to the natural frequencies of the tall buildings 
existing in this area, suggests a local selective 
amplifying effect [21]. The existence of softer 
layers at Caraballeda may be responsable for 
the more severe damage that occured in this area. 
However this important problem deserves further 
study. 


4.3 — Behaviour of buildings 


Fig. 14 shows that at Palos Grandes and Alta- 
mira area the buildings that collapsed or suffered 
important damages had coefficients of resistance 
below about O 08. 

The fact that Petunia II building (Table I, ref. 
no. 8), having a coefticient of resistance of 0.10, 
underwent important damage can in part be 
explained by its cantilever behaviour which 
corresponds to a lower ductility than the one 
corresponding to frame behaviour. 

In the Caraballeda area the buildings that 
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collapsed or suffered important damages had 
coefficients of resistance below about 0.22. In 
Naiguatá and Playa Grande this limit decreases 
to about 0.12. 

The influence of the natural frequency on the 
amount of damage is not directly apparent in 
fig. 14. This may be due to the relatively small 
range óf natural frequencies in the structures 
and also to the fact that the influence of this 
factor is not as important as indicated by the 
proportionality between seismic coefficient and 
natural frequency. 

Fig. 15 relates the measured natural periods, 
T, in longitudinal and transverse directions to 
the mumber of stories, N. The mean value of 
variation, V, being 35º/,. 


Correlations of the type T(s) = K 7 , where 


K is a constant, H the height and B the depth 
(m) of the buildings were also considered the 
following results being obtained. 


Number ” 

HjB of cases x V (to) 
= 7) | 29 0.082 14 
< 2 19 0.123 | 22 
all 48 0.098 29 


Particular reference is due to the very good 
behaviour of Plaza 1 building (Table IV, ref. 
no. 42) having a reinforced concrete shear walls 
structure. 

A description of some accidents due to faulty 
details is not included in the present paper but 
can be found in several of the indicated refe- 
rences. 


5 — CONCLUSIONS 


The main conclusions derived are the following. 


5.1 — The intensity of the earthquake was very 
different in the various areas of Caracas and 
nearby locations. The following approximate 
estimates of the power spectral density of 
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. . 2 — 4 
acceleration expressed in cm'S Hz, can be 
indicated : 


Caraballeda . . . ... 2 1500 
Playa Grande — Naiguat T 700 
Palos Grandes — Altamira . .7 500 
Cajigal = 750 


5.2 — The variation of seismic intensity seems 
mainly due to different soil conditions. Even so 
present information does not permit a quantita- 
tive assessment of the relationship between the 
observed intensity and the soil conditions. 


5.3 — The comparison between the computed 
coefficients of resistance and the observed 
damages shows that the buildings that collapsed 
or suffered important damage had coefficients of 
resistance below about 0.22 in Caraballeda, 0.12 
in Naiguatá and Playa Grande, and 0.08 in Palos 
Grandes — Altamira. 


54—In the observed cases, the influence of 
the natural frequency of vibration on the relation 
between coefficients of resistance and amount of 
damage does not appear to be important. 


5.5 — Due to the fact that most structures 
were formed by frames, the influence of the type 
structure and of the correspondent ductility on 
the relation between coefficients of resistance 
and amount of damage was not clearly detected. 
It would be of great value to obtain further data 
on this matter. 


5.6 — The collapses can in general be explained 
by the small values of the coefficients of resis- 
tance and are not attributed to poor construction 
or faulty detailing. The coefficient of resistance 
are particularly low for those structures that in 
one direction had no beams but only slender 


flat slabs. 


5.7 — The presented results may contribute in 
assessing the adequacy of seismic factors usually 
adopted. Attention must be paid to the fact that 
the indicated coefficients of resistance can only 
be directly compared with the seismic factors 
indicated by those codes that consider ultimate 
design and load factors equal to 1. 
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NOTA SOBRE A INFLUÊNCIA NAS CONSTRUÇÕES DAS 
VIBRAÇÕES DOS TERRENOS PROVOCADAS 
POR EXPLOSÕES” 


RESUMO 


Neste trabalho apresenta-se um estudo da bibliografia 
mais recente sobre a matéria em epígrafe. 

Indicam-se os critérios de segurança que têm sido 
aplicados nos E.U.A., Canadá e Suécia concluindo-se que 
o critério baseado na velocidade de vibração das particulas, 
hoje o mais adoptado, é uma medida da deformação 
experimentada por uma construção fundada num terreno 
onde se propaga uma onda de vibração, 

O limite geralmente admitido para a velocidade máxima 
de vibração das partículas dum terreno é de 50 mm/s 
(2 in/sec), abaixo do qual é de prever que não haverá 
danos na construção. 

Apresentam-se também alguns resultados obtidos no 
L. N. E. €., concluindo-se ser necessário estudar para cada 
caso, o modo de propagação da vibração provocada por 
explosões, 

Finalmente, com o objectivo de se iniciar um trabalho 
com boa técnica, apresentam-se alguns elementos que permitem 
numa primeira aproximação calcular a carga de explosivo 
a empregar, 


1— INTRODUÇÃO 


Pode dizer-se que desde as primeiras aplicações 
de explosivos em trabalhos mineiros e de enge- 
nharia civil, se levantou o problema dos efeitos 
nas edificações das vibrações dos terrenos e do 
ar originadas pela sua detonação. Este problema 
é especialmente importante, quando se reali- 
zam desmontes de rocha empregando explosivos, 
nas proximidades de zonas habitadas. 

Não existe ainda, porém, um critério de segu- 
rança igualmente seguido na generalidade dos 
países, definidor do máximo nível de vibrações 
provocadas por explosões que não deve ser ultra- 
passado sob perigo de ocorrerem danos nas edifi- 
cações situadas nas suas proximidades. 

Este problema tem sido objecto de vários estu- 
dos especialmente na Suécia, E.U.A. e Canadá, 


por JOAQUIM MOURA ESTEVES** 


SYNOPSIS 


This paper is a study of the most recent literature on 
the subject in reference, 

The safety criteria in use in the US.A, Canada and 
Sweden are indicated, it being concluded that the criterion 
based on the velocity of propagation of particles, the most 
used today, is a measurement of the strains in à structure 
founded on soil where the vibration wave progresses. 

The usual limit velocity of the particles in the soil is 
50 mm/s (2 in/sec), below which no damages are to be 
anticipated in the structures, 

Some results obtained atthe L.N.E.C. are also presented, 
it being concluded that the propagation characteristics of 
blast vibrations should be studied in each particular case, 

Finally in order to facilitate the application of safe 
methods, some data are presented by means of which the 
amount of explosive to be used can be calculated in a firet 
approximations, 


subsidiados por entidades estatais, companhias 
de seguros, fabricantes de explosivos, cons- 
trutores, etc. Entre nós, que se saiba, não 
tem havido estudos específicos do problema e 
o que se tem feito tem tido apenas carácter 
qualificativo. 

Sucintamente referem-se a seguir os critérios 
de segurança aplicados a edificações, desde 1942, 
nos E.U.A., Canadá e Suécia, países que têm 
estudado, como se disse, mais profundamente a 
influência mas construções das vibrações dos 
terrenos provocadas por explosões. Aquela data 
marca o início das investigações realizadas no 
Bureau of Mines dos E.U.A. com vista ao esta- 
belecimento dum critério de segurança definidor 
do máximo nível de vibrações originadas por 
explosões, que as edificações podem suportar sem 
nelas se verificar dano. 


* 2.95 Jornadas de Engenharia e Arquitectura do Ultramar — Luanda, Maio e Junho de 1969, 
** Engenheiro Civil, Especialista do Laboratório Nacional de Engenharia Civil. 
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Neste trabalho analisa-se o critério de segu- 
rança que, com base naquelas investigações, 
parece dever seguir-se entre nós, quando não se 
realizam estudos específicos, suficientemente por- 
menorizados, e apontam-se regras práticas para 
o cálculo da carga de explosivo a utilizar em cada 
detonação e para certos tipos de terreno. 


2 — CRITÉRIO DE SEGURANÇA 


O primeiro relatório das investigações efectua- 
das pelo Bureau of Mines, em 1942 [1], referia 
o critério de segurança relativamente à rotura ou 
dano duma construção baseado nas acelerações 
experimentadas pela estrutura, de modo que se se 
atingisse um nível de vibração inferior a 0,1 g 
(um décimo da aceleração da gravidade) isso 
significava não existir perigo para a edificação; 
mas para níveis de vibração correspondentes a 
valores de 0,1 g a 1 g haveria que tomar pre- 
cauções; acima de 1 g a situação era de perigo. 

Pouco depois desta data, em 1949, F. J. Crandel 
[2] estabelece um critério de segurança, duma 
construção, especificamente baseado no nível de 
vibração alcançado pelo terreno onde se situava 
a construção. Este critério de Crandel é expresso 
em termos dum coeficiente energético, definido 
pelo cociente entre o quadrado da aceleração e o 
quadrado da frequência da vibração ; a aceleração 
é expressa em ft/s* e a frequência em c/s. Para os 
coeficientes energéticos inferiores a 3, a situação 
era considerada de segurança; entre os valores 
3 e 6 teria de se tomar alguma precaução ; acima 
do valor 6 a situação era de perigo. 

Depois, até 1960, observou-se um crescente 
interesse por estes problemas, mas, no entanto, 
não se encontrou um critério de segurança geral- 
mente aceite. Por exemplo, nos E.U.A., adopta- 
vam-se os mais variados critérios. Assim, no 
estado de Massachussetts, considerava-se o valor 
1 do critério energético como limite de segurança, 
enquanto no estado da Pensilvânia se considerava 
o valor de 0,03 polegadas (critério de desloca- 
mento) para esse limite de segurança. 

Noutros países, especialmente na Suécia e no 
Canadá, começou a adoptar-se um critério de 
segurança baseado na velocidade de vibração das 
partículas do terreno onde se situava a constru- 
ção e escalonado do seguinte modo: na Suécia, 
em 1957 [3], para velocidades inferiores a 
70 mm/s (2,8 in/sec) não se observavam danos; 
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deste nível de velocidade ate 110 mm/s (4,3 in/sec) 
poderiam surgir fendas no estuque e mesmo que- 
das deste; até 160 mm/s (6,3 in/sec) surgiriam 
algumas fendas na estrutura; para 230 mm/s 
(9,1 in/sec) seriam já graves as fendas que 
surgissem. No Canadá, em 1959 /4 |, foi proposto 
um critério baseado também na velocidade de 
vibração das partículas do terreno, correspon- 
dendo o nível de 50mm/s (2 in/sec) ao limite 
considerado seguro; entre 50 e 100 mm/s (2 e 4 
in/sec) haveria que tomar precauções; acima 
de 100 mm/s (4 in sec) ocorreriam grandes danos 
ou rotura. 

Em face deste panorama, de grande contro- 
vérsia, sobre qual o critério de segurança que 
deveria ser adoptado, e tendo em consideração 
o incremento do número de problemas criados 
pelo desenvolvimento do emprego de explosivos 
e atendendo ainda ao desenvolvimento verificado 
na tecnologia do desmonte a fogo, decidiu 
o Bureau of Mines (E.U.A.) estudar de novo 
o problema das vibrações originadas por este 
tipo de desmonte e os seus efeitos nas edifi- 
cações. 

O mais importante desses estudos, feito com 
aplicação de métodos estatísticos [5], teve por 
base, especialmente, os trabalhos realizados pelo 
próprio Bureau of Mines, pelos investigadores 
suecos Langefors, Kihlstrôm e Westerberg, e pelos 
investigadores canadianos Edwards e Northwood 
(1, 3 e 4], pois, embora tivessem sido estudados 
47 trabalhos, apenas os realizados por esses 
investigadores continham valores de amplitudes 
e frequências de vibrações, originadas por 
explosões, relacionados com os danos por elas 
provocados em edificações. Desse estudo consta- 
tou-se, em todos os 124 casos analisados, que 
era a grandeza velocidade de vibração das 
partículas do terreno que melhor caracterizava 
a interdependência entre a vibração do terreno 
e os danos por ela provocados nas edificações 
sobre ele construídas. Esta constatação é bastante 
significativa uma vez que os dados que a origi- 
naram foram obtidos por investigadores traba- 
lhando independentemente uns dos outros e utili- 
zando variados aparelhos de medida em locais 
diversos, onde existiam diferentes edificações. 
Parece por isso aplicável, mesmo para diversas 
condições físicas, o critério de segurança contra 
danos ou rotura de edificações baseado na velo- 
cidade de vibração das partículas do terreno. 
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Como o estudo se baseou em métodos estatís- 
ticos apenas foi possível definir uma fronteira 
correspondente ao nível de vibração de 50 mm/s 
(2 in'sec), abaixo do qual não é de esperar que 
ocorra dano na edificação. Este limite não se 
afasta muito dos limites estabelecidos por outros 
investigadores. Por exemplo, o coeficiente ener- 
gético 3, de Crandell, corresponde a uma velocidade 
de 84 mm/s (3,3 in/sec); para Langefors, 
Kihlstrôm e Westerberg esse limite é de 70 mm/s 
(2,8 in/sec); e para Edwards e Northwood é de 
50 mm/s (2 in/sec). 

Aquele estudo mostra ainda que o nível de 
190 mm/s (7,6 in/sec) corresponde a 50º/o da 
probabilidade de ocorrerem grandes danos (queda 
de estuque, graves fendas) e que a probabilidade 
de ocorrerem pequenos danos (ocorrência de finas 
fendas no estuque, abertura de fendas já existen- 
tes) é ligeiramente inferior a 50º/ para uma 
velocidade das partículas de 137 mm/s (5,4 in/sec). 

O critério de segurança baseado na velocidade 
de vibração das partículas, que é hoje o adoptado 
duma forma geral, parece ser razoável, pois, como 
observam Ambraseys e Hendron [6], a máxima 
velocidade radial Vr duma partícula é propor- 
cional à máxima deformação unitária radial 
conforme mostra a equação (1): 


r 


== — (1) 


em que c é a velocidade de propagação da vibra- 
ção no terreno. 

Assim, a velocidade da partícula é uma medida 
da deformação do terreno, sendo por isso de 
esperar que tal velocidade seja uma medida da 
deformação e dano observado numa construção 
fundada nesse terreno, onde se propaga uma 
onda de vibração. 

Considerando então aceitável, em face do conhe- 
cimento actual do problema, que é a grandeza 
velocidade de vibração das partículas do terreno 
que melhor caracteriza a interdependência entre 
a vibração do terreno e os danos que eventual- 
mente surgem nas estruturas sobre ele edificadas 
e que não se deve ultrapassar o nível de veloci- 
dade de 50mm/s (2 in'sec), haverá agora que 
saber qual a máxima carga de explosivos que 
pode ser detonada sem que, portanto, se originem 
danos. 
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3— CARGA DE EXPLOSIVOS A UTILIZAR 


A carga máxima de explosivos a empregar 
depende não apenas do nível de vibração local 
que origina a sua detonação mas também da 
distância da explosão ao edifício e das proprie- 
dades dinâmicas do meio onde se transmite a 
vibração. 

Com base nas medições de velocidades máximas 
de vibração das partículas feitas em 19 locais de 
pedreiras, W. I. Duvall, €. F. Johnson. A. V. €. 
Meyer e J. F. Devine [7] estabeleceram a seguinte 
equação da vibração : 


= KQ'D”" (2) 
em que: 


V — velocidade máxima de vibração da parti- 
cula 
K — constante dependente do meio de propa- 
gação 
Q — carga de explosivo por rebentamento 
D — distância entre os pontos de explosão e 
de observação 
b e n— expoentes que variam com a compo- 
nente da velocidade (radial, vertical, 
transversal) 


Ambraseys e Hendron [6] concluem da aná- 
lise dimensional que a velocidade é função da 
relação D/Q!" entre a distância D e a raiz cúbica 
da energia ou carga explosiva Q. Embora 
alguns autores como Devine [8] apresentem 
valores de velocidade em função da relação 
D/0"*, parece ser, efectivamente, de recomendar 
a relação D Q!” visto que sendo baseada na aná- 
lise dimensional permitirá extrapolar a regra ou 
lei estabelecida a partir dos elementos existentes 
ou obtidos localmente. 

Seguindo de perto Ambraseys e Hendron [6] 
apresenta-se a figura 1, onde se observa a variação 
dos valores da velocidade máxima da partícula 
em função dos valores da relação D/Q!”. Para 
valores superiores a 10ftI/b'” recomenda-se 
a equação (3) 


Vr == 360 in/sec (D/Q'”) gi (3) 


sendo Vr a velocidade radial, com D expressa 
em pés e Q em libras. Para menores valores da 
relação D/Q! recomenda-se a equação (4) 


Vr = 6000 in/sec (D/0!")— 2º (4) 
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correspondendo esta lei a uma maior atenuação 
da velocidade máxima do que para o caso anterior, 
devido a comportamento inelástico do terreno. 
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radiais, das partículas do revestimento do túnel 
do Rossio, obtidos pelo L. N. E. C. aquando dos 
trabalhos de abertura, com emprego de explosivos, 
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Fig. 1 — Valores de campo de velocidades máximas das partículas em função da 
relação D/O! (origem: Ambraseys and Hendron [6]) 


Nesta mesma figura estão implantados alguns 
pontos correspondentes a valores de velocidades 
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da Avenida Gulbenkian sobre aquele túnel [9]. 
Embora estes resultados sejam representativos 
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apenas dum dado local, verifica-se que os corres- 
pondentes pontos se dispõem segundo uma lei 
diferente das acima referidas, o que revela que 
elas não devem ser indistintamente aplicadas. 

Deverá, pois, efectuar-se para cada caso os 
estudos do modo de propagação da vibração 
provocada por explosões, procurando-se estabe- 
lecer uma relação entre a velocidade máxima de 
vibração das partículas do terreno junto à edifi- 
cação e o cociente D/Q!* entre a distância D do 
ponto de detonação ao ponto de recepção da 
onda vibratória e a raiz cúbica da carga de 
explosivo. 
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ds TANCE 
Fig. 2 — Carga de explosivo O em função da distância 
D para diversos valores da relação 0/DI!5 
(origem: U. Langefors and B. Kihlstrôm (ro]). 


No entanto, e a título de início dos trabalhos, 
recomenda-se a utilização do ábaco da figura 2 
e do quadro | estabelecidos por Langefors e 
Kihlstrôm para se calcular a carga de explosivos 
a empregar [10]. 

Estes mesmos autores apresentam também 
num outro quadro (quadro II) o tipo de dano 
ou rotura que pode ocorrer em habitações tipo 
moradia, fundadas em terrenos de natureza 
diversa, relacionando-o com o nível de veloci- 
dade de vibração das partículas desses terrenos. 
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Daqui ressalta a importância do conhecimento 
da velocidade de propagação das ondas elásticas 
em cada tipo de terreno onde se processa a 
utilização de explosivos. 

Efectivamente, o valor de 50mm/s (2 in/sec), 
estabelecido como limite de segurança para a 
velocidade de vibração das partículas do terreno, 
poderá ser elevado até 70 mm/s (2,8 in/sec) se 
a construção está fundada num terreno rochoso, 
mas terá, eventualmente, que ser reduzido se 
a edificação se encontrar implantada num terreno 
com baixa velocidade de propagação das ondas 


; Vr 
elásticas. A expressão (1) :r = mostra que 


ao pretender-se manter o valor da extensão Er 
da partícula do terreno, e portanto a sua defor- 
mação aquando da utilização de explosivos, a 
menores valores da velocidade c de propagação 
das ondas elásticas correspondem menores valo- 
res para a velocidade Vr de vibração das par- 
tículas, e vice-versa. 

Na figura 3 apresenta-se um exemplo de 
registo obtido pelo L.N.E.C. no trabalho já 
referido do controle do nível de vibrações pro- 
vocadas no revestimento do túnel do Rossio pela 
abertura a fogo da Avenida Calouste Gulbenkian. 
Foram registadas velocidades de vibração das 
partículas em três direcções ortogonais cujas 
grandezas foram calculadas com base em expres- 
sões previamente determinadas em ensaios lato- 
ratoriais de calibração da aparelhagem de medida. 
As explosões foram efectuadas com cápsulas 
retardadas de 0,5 milisegundo de modo a evitar- 
-se a acumulação de efeitos. Facilmente se observa 
no registo que não se mantém aquele intervalo 
de tempo entre as diferentes explosões. Isto é 
devido não só a imprecisão de fabrico das 
cápsulas microrretardadas, mas também ao dife- 
rente modo de propagação da onda vibratória de 
ponto para ponto de explosão. 

Finalmente e como se depreende do exposto, 
o estudo da influência das vibrações dos terrenos 
provocadas por explosões em edificações neles 
implantadas constitui vasto campo de investi- 
gação em diversos países especialmente nos 
E.U. A. Dada a importância do problema julga- 
-Se conveniente que se promova entre nós tam- 
bém o seu estudo, pelo menos até ao ponto de 
se poder quantificar o nível de vibrações origina- 
das pela utilização de explosivos nos centros 
urbanos. 
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Fig. 3— Exemplo dum registo de várias explosões 


QUADRO 


Quantidades de explosivo para diferentes valores da relação Q/D a 


Distâncias Q tg) 

(m) op oo 0,06 0,12 01,25 0,50 1,0 
0,4 0,07 0,14! 0,28 
1,0 0,03 0,06 0 0,24 0,5 1,0 
14 (1,05 0,10 0,2 0,1 0,8 1,6 
2,0 0,09 0,2 0,4 0,7 1,4 2,8 
28 0,15 (1,8 0,6 1,2 2,4 4,8 
4,0 0,25 0,5 1,0 2,0 4,0 8,0 
6,0 0,5 0,9 1,8 3,1 7,5 | 15 

10,0 1,0 2 4 S 16 13 

14 ti 3,4 7 14 28 56 

2) 28 5,8 1! 29 44 SS 

28 5 10 20) 40 80 160 

40 8,5 17 35 70 140 280 

60 15 30 60 I20 249 480 

100 32 64 130 2"0) 20 | 

140 50 100 210 429 810) 

200 90 | 180 360 720 1410 

280 150 300 600 1200 2400 

400 250 500 1,000 | 2000 4000 

600 500 900 1,800 | 3700 

1000 1000 2000 4,000 | 8000 
1400 1700 3400 | TO 

2000 2800 | 5600 | 11,000 

2800 5000 | 

4000 | 8500 | 


Origem ; U. Langefors and B. Kihlstrôm [10]. 
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QUADRO Ii 


Tipo de dano em residências género moradia implantadas em diferentes terrenos 


Natureza Areia, cascalheira 
do argila abaixo do 
terreno nivel freático 


| Moreia, xisto, 
calcário brando 


Calcário duro, grés 
quartzítico gnaisse 
granito, diabase 


Tipo de dano 


E — as | EO À LE 


Velocidade de propagação 


c (m/s) 1000-1500 2000-3000 
Velocidade 18 35 
de Bil do 
vibração V 40 80 
(mms) 60 115 


4500-6000 
70 Não se observam fendas 
no Fendas insignificantes 
160 Fendilhação 
230 Grandes fendas 


Origem: u. Langefors and B. Kihlstrôm [10]. 
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ELEMENTOS SOBRE A PRODUÇÃO E O CONSUMO DE ENERGIA 
NA REDE ELÉCTRICA NACIONAL *) 


— Breve nota mensal 


Em Junho registaram-se ainda precipitações, que 
embora modestas em valor absoluto, se situaram em 
valores da ordem das correpondentes médias. 

Em consequência das elevadas precipitações ocor- 
ridas nos últimos dias de Maio, que se reflectiram nas 
afluências «internas» do mês de Junho, e do caudal 
bastante favorável transitado em Miranda, as afluências 
ao conjunto do sistema produtor, durante o mês de 
Junho, atingiram cerca de 571 GWh valor este que tem 
uma probabilidade de ser excedido da ordem dos 8º/, 
e a que corresponde um coeficiente de produtibilidade 
hidroeléctrica de 1,67. 


II — Elementos gerais (GWh) 


a) Mensais Variação 
1968 | 1969 | “Vo 
Produção hidráulica (Ph) . « | 589,2 | 471,9 + 21 
Produção térmica (P+)..... .| 838 0,0 !—100 
Produção total (PT)......... 427,5 | 471,5 + 10 
Energia recebida de empresas 
não pertencentes ao RNU (Er) | 1,4 2,7 + 93 
Exportações (Ex).,....,.ce.. 3,3 2,8/— 15 
Importações (1)... .-cc.... ass 1,5 2,2 |-+- 47 
Saldo importador (Sp)....... —1,S | —0,6/+ 67 
Cons. em bomb, hidroelé. Cb) (1) 0,0 0,U 0 
Produção para con- (1) 
sumos perman. (Pep)..... 878,1 |, 419,4 /4-10,9 
Produção para con- 
sumos não perman. (Penp) ..., 490, 54,2 +11 
O OGRO 4271 | 473,6|+ 11 
Coeficiente de hidraulicidade 0,74] 1,67) — 


NOTA: 


(1) O aumento percentual da produção para consumos permanen- 
tes, tendo em conta a incidência dos domingos e dias especiais 
e respectivamente de 11,9 e 9,40» 


III — Diagramas de carga dos dias característicos 


4.º feira; 


—tem —— 


1968 | 1969 
Produção hidráulica  (P ) MWh [13121 17821 
Produção térmica (Po) MWh | 2290 0) 
Produção total (Pr) MWh | 15411 17821 
Trocas com f Export.  (|x)MWh 75 130 

Espanha À Import. (1 )MWh OU VU 
Cons. em bomb. hidroel. (Cb) MWh UV | O) 
Prod, p.* cons. perm. (Pep)MWb|14102 | 16024 
Prod.p.*cons.nãoperm.(Penp)MWb| 1284 17957 


TOTAL  Pr+(Ex)MW5 [14586 (17761 - 
o Potênciamax. MW| Sol As. 
“gps Ex)! Potência min, MW 393 46 
s e Pr+t ne Utiliz. da ponta horas 18,0 18,2 
E e “| Factor de carga 0,70 0,76 
pre Potência máx. MW TOY 901 
Es p Potência min. MW 338 394 
dis e Utiliz. da ponta horas 17,6 17,8 
Factor de carga 0,73 0,74 
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Iv — Energia armazenada nas principais albufeiras 


No fim do mês 


Albufeiras: 


GWh 91 43) 
Alto Rabagão +... ... «| 969,9 99,6 
Paradela. . cc cv 00. 218,7 98,5 
Venda Nova +. «cc» «o| 1199 93,7 
Salamonde . « «cc cv... 27,3 98,9 
Caniçada ..ccrcvoo 10,4 31,4 
VEMOS cs ds» asa md 99,0 98,4 
Cabril. ucse coco «o +| 'BB85 99,8 
Castelo do Bode . . . «« «| 1541 34,9 
Guilhofrei . . . cc cv... 6,6 79,5 
Lagoa Comprida .« +... «-«|(02B7Tb 98,0 
Rasta Eneida usa mes 58,0 94,2 
Pracanã » cc creo vo. 10,7 82,9 
Póvoã » cc oc voo» «| (10,6 76,0 
À com A, Rabagão . . .| 2050,9 91,1 
Total sem A, Rabagão. . 1091,4 95,0 


NOTAS 
(1) Coeficiente de enchimento em energia das albufeiras, 


(2) Inclui 3,0 GWh armazenados em Vale do Rossim no início do 
mês e 3,3 GWh no fim do mês, 


(3) Inclui 2,7 GWh armazenados no açude do Poio no início do mês 
e 2,7 GWh no fim do mês. 


(*) Elementos extraídos das estatísticas mensais do Repar- 
tidor Nacional de Cargas (R. N, C.). As produções e 
os consumos das empresas do R.N.C. representam 
cerca de 94º/y dos totais do País. 
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CONDUÇÃO DO CALOR EM REGIME NÃO PERMANENTE, 
DE UM FILME PARA DOIS MEIOS SEMI-INFINITOS 


RESUMO 


Estuda-se a evolução no tempo, da temperatura de 
um filme condutor, quando aquecido por efeito Joule, 
no caso do calor gerado ser transmitido apenas por 
condução aos dois meios semi-infinitos que o envolvem. 

Admite-se como constante a quantidade de calor 
gerado por unidade de tempo, em cada unidade de 
superficie, e considera-se a existência da inércia térmica 
do filme. 


0. — NOMENCLATURA 


por SILVESTRE DIAS ANTUNES 


Bolseiro da Universidade de Luanda 
no Núcleo de Estudos de Eng. Mecânica * 


SYNOPSIS 


It is developed a theoretical method of calculation of 
the temperature from a film heated by an electrical cur- 
rent, with uniform power Q per unit of surface. 

The heat is transferred from the film to the two semi- 
Infinite mediums by conduction, considering the ther- 
mal capacity of the film. 


A — capacidade calorífica do filme, referida à unidade de superfície (inércia térmica do filme) 


B — constante definida em (4.12) 
C, (ns: inteiros) — constantes de integração 
2 


ax 
erÊ x = J exp (—u?) du — «error function» 
Va do 


erfcx=-1 — erf x 
exp ( ) — função exponencial 
£[ ]—transformada de Laplace 


£2"![ |] —transformada de Laplace inversa 


p — variável complexa da transformada de Laplace 


Q — quantidade de calor gerado no filme, em cada unidade de superfície e por unidade de tempo 


t —tempo 

T —temperatura 

« — coeficiente de difusidade térmica 

k — coeficiente de condutibilidade térmica 
à — espessura do filme 


O (x, t) =T— To — acréscimo de temperatura desde o instante inicial, até ao instante t 
O (x,p)= [0 (x,t)] — transformada de Laplace de O (x,t) 
9º (x, t)=lim O (x,t) — acréscimo de temperatura desde o instante inicial até ao instante t, quando 


A -»0 


se despreza a existência de capacidade calorífica (A + 0) 


índices: 


o — instante inicial 
1, 2, 3 —- meios que constituem o sistema 


* Laboratório de Termodinâmica Aplicada — Departamento de Eng. Mecânica — I. S. T. 
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1. — INTRODUÇÃO 


No âmbito do programa de actividades do Laboratório do N.E.E.M., no estudo de técnicas 
de medida da condutibilidade térmica de materiais, surgiu o problema da transmissão de calor, por 
condução, em regime transiente, de um filme para dois meios semi-infinitos. Como, na bibliografia 
consultada se não encontrou qualquer solução para o problema estudado, apresentam-se a seguir as 
conclusões a que se chegou. 


2. — MODELO. HIPÓTESES 


Um filme condutor, plano e infinito, separa dois meios semi-infinitos, homogéneos e isótropos, 
de materiais de natureza diferente. 

Pretende estudar-se a condução do calor gerado no filme por efeito Joule, para os meios que 
o cercam, para o que se considera um modelo unidimensional, de direcção característica perpendicular 
ao plano do filme. 

O sistema constituído pelo filme e pelos dois meios semi-infinitos é fisicamente caracterizado 
pelas seguintes hipóteses: 


H.1 — no instante inicial a temperatura é uniforme em todo o sistema; 

H.2 — os coeficientes térmicos dos materiais são independentes da temperatura; 

H.3 — o calor gerado no filme é transmitido apenas por condução aos dois meios que o envolvem, 
nos quais se supõe não haver fontes de calor ; 

H.4 — nas superfícies de contacto não existe resistência térmica (contacto térmico perfeito); 

H.5 — dada a pequena espessura do filme, e a grande condutibilidade térmica do material que 
o constitui, relativamente à dos materiais que o envolvem, admite-se como uniforme, 
em cada instante, a temperatura em todo o filme (hipótese esta que equivale a admitir 
como infinita a condutibilidade térmica do material do filme) ; e, 

H.6 — a quantidade de calor gerado por efeito Joule, por unidade de tempo e de superfície, 
é constante, a partir do instante inicial. 


3. —- FORMULAÇÃO MATEMÁTICA 


No sistema considerado a evolucão da temperatura em cada instante e em cada ponto rege-se 


pela equação do calor : 
d e) 
o = + Ip 6 4 O (3.1) 
dt pc GC 

em que: 


O) =T — To — acréscimo de temperatura do ponto considerado desde o instante inicial até ao 
instante t 
k — coeficiente de condutibilidade térmica 
o — massa específica 
c — calor específico 
k é iuá di à , 
a — — — coeficiente de difusividade térmica 
pc 
Q — calor gerado por unidade de volume e de tempo. 


Considere-se o sistema de coordenadas definido na figura 1, então (3.1) escrever-se-á: 


2 
elCui ki dd Cd) (3.2) 
dt dx 
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Pretende determinar-se 9 (x,t) condicionada pelas seguintes condições : 
CONDIÇÕES INICIAIS 

C1— (O (x,0) = o qualquer que seja x (hipótese H.1) (3.3) 
CONDIÇÕES AOS LIMITES (razões de ordem física) 


C.2— O (x,t) é contínua em todo o domínio de x, e em particular para x=: 


lim O (xt) = Hi, O (xt) = O (os) (3.4) 


R-+ 0 X —> O 


C3- O (xt) é finita em todo o domínio de x, qualquer que seja t F oo, e em particular 
para x = + 00: 


lim O (xt) é finito se t E co (3.5) 
lim O (xt) é finito se t & oo (3.6) 
x —» +00 


Meio (2) 


meio (1) meio (2) 


a, K, a, K;, 


Fig. 1 


ER Ta 


Fig. 2 
Consideremos o elemento de volume representado na figura 2, constituído por um cilindro 


de geratrizes perpendiculares ao plano do filme e de secção recta S. Integrando a equação (3.1) no 
volume V definido pelo cilindro considerado, obtém-se a forma integral da equação do calor: 


d (3) -. . 
f oc — dV a k lap BdV+ | Q dV (3.7) 
vo dt V dv 


e sea altura do cilindro se for reduzindo até que este coincida com um sector do filme, vejamos a 
forma que tomam cada um dos termos de (3.7). 


Da definição de inércia térmica, vem: 


lim c—dV= as(S) 3.8 
ah.,8 vo o dt = 0 
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